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ABSTRACT 

Dark Matter supposedly dominates the extragalactic, yet no totally dark structure 
of ga lactic proportions has ever been convincingly identified. Earlier ( Min chin et alJ 
2005) we suggested that VIRGOHI 21, a 21-c m source we found in the Virgo Cluster 
at Jodrell Bank using single-dish observations l|Davies et a l. 2004), was probably such 
a dark galaxy because of its broad line-width (~ 200 kms -1 ) unaccompanied by any 
visible gravitational source to account for it. Now we have managed to image VIR- 
GOHI 21 in the neutral-hydrogen line, and indeed we find what appears to be a dark, 
edge-on, spinning disc with the mass and diameter of a typical spiral galaxy. Moreover 
the disc has unquestionably interacted with NGC 4254, a luminous spiral with an odd 
one-armed morphology, but lacking the mas sive interactor invariably responsible for 
such a feature. Published numerical models ijVollmer. Huchtmeier fc van Driell f2005) 
of NGC 4254 call for a close interaction ~ 10 8 years ago with a perturber of ~ 10 11 
solar masses. This we take as completely independent evidence for the massive nature 
of VIRGOHI 21. 
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1 INTRODUCTION 

Cold Dark Matter (CDM) models of galaxy formation 
predict many mo re dark matter halos than are ob- 
served as galaxies feauffmann. White fc Guiderdonilll993t 
iMoore et al1ll999l : iD'Onghia fc Lakdl2004lh Whether stars 
can form in dark matter halos depends critically on the frac- 
tion of baryons (mj) that can be trapped and retained by 
each halo (for previous discus sions about the existence and 
formation of dark galaxies see ( Ji menez et al Jl997tlHaw]dnl 
ll997tlVerde et alJl20oA l. Some of those that retain a small 
fraction of their original baryons (vtid < 0.05) are able to 
form stable gaseus disks, but because of the low gas den- 
sities they are not able to form stars. The physical state 
of the gas in the disc is very dependent on its density and 
temperature with the latter influenced by the intensity of 
any ionising background. Models predict that galaxies can 
form with gas column densities that prohibit star formation 
yet provide some self shielding from the ionizing background 
(Davies et al. in prep.). Such dark galaxies are potentially 
detectable by blind 21cm surveys of the sky. 

Objects detected at 21cm but with no optical coun- 



terparts have been known about for many years; these in - 
clude high ve locity clouds i jWakke r fc van Woerdenlll997f) . 
the Leo Ring JSchneider et alll983ft. an d various gas clouds 
close to bright galaxies (iKilborn et alJ l200Ct iBovce et alJ 
l200lt iRvder et alJl200llh However, none of these objects 
have the characteristics of a galaxy, i.e. detectable emission 
over galaxy sized spatial scales and a velocity structure con- 
sistent with a rotating and gravitationally bound disc. 



In a previous paper (iMinchin et alJl2005h we described 
Jodrell Bank and Arecibo observations of VIRGOHI 21 an 
HI source discover ed during a survey of the Virgo cluster 
jDavies et al.ll2004ft . VIRGOHI 21 has all of the HI charac- 
teristics of a typical rotating disc galaxy, but no detectable 
optical emission. Our conclusion was that VIRGOHI 21 was 
an extremely promising candidate for the first dark galaxy. 
In this paper, we present high resolution HI observations 
of this source, which add to the conclusions drawn in the 
earlier paper that this is, indeed, a massive dark galaxy. 
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2 OBSERVATIONS AND ANALYSIS 

The new data were taken in March 2005 at the Westerbork 
Synthesis Radio Telescope (WSRT) in two full 12-hour syn- 
theses and reduced using the Miriad package. The data 
were flagged for shadowing and on two of the fourteen 25- 
metre antennae one polarisation was flagged due to problems 
with the gain. A spectral bandwidth of 10 MHz covered the 
velocity range 930 - 3070 kms" 1 . Removal of the noisier end 
channels left 230 useful channels of width 8.2 kms -1 each, 
giving a velocity resolution of 10 kms -1 over the range 980 
to 2890 kms . Continuum removal was carried out in the 
UV plane using Uvlin. The standard source 3C147 was used 
for calibration. Cl eaning used a robust setting of 1, close to 
normal weighting feriggsll995ft . The cleaned cube was gaus- 
sian smoothed spatially and Hanning smoothed in velocity, 
and then used for a second-pass deeper cleaning which gave 
the final cube used in the analysis. The synthesised beam 
was 99" x 30" in size (extended North-South) and the noise 
was 0.3 mjybeam -1 channel -1 , giving a 5<r column-density 
to sources 25 kms -1 wide of 2 x 10 19 Hydrogen atoms cm -2 . 

Fig. la shows a neutral Hydrogen (H i) contour map of 
the field superimposed on a negative optical image. VIR- 
GOHI 21 is the elongated structure in the centre (which is 
at about 2000 kms -1 ). A faint bridge can be seen stretch- 
ing down to the prominent spiral NGC 4254 (2400 kms -1 ) 
while the other two sources, NGC 4262 (1500 kms -1 , upper 
left) and the faint galaxy 'C (1750 kms -1 , immediately to 
the left of VIRGOHI 21) appear unconnected. Fig. lb shows 
the velocity-declination projection of the full 3 dimensional 
data cube. Now NGC 4254 is at the bottom right while VIR- 
GOHI 21 is the angular structure in the centre. C is to its 
left. Far more detail can be seen on an animation which is 
available at http://www.astro.cf.ac.uk/groups/galaxies for 
instance the bridge is clear and obvious, as is the lack of any 
connection between VIRGOHI 21 and either NGC 4262 or 
C. The apparent alignment of NGC 4262 with VIRGOHI 21 
and NGC 4254 in Fig. lb is a consequence of the particular 
projection shown 

The ~ 25 arcmin filamentary bridge stretches from the 
low-velocity (western) edge of NGC 4254, falling gently in 
radial velocity from 2250 kms -1 at declination +14°20' to- 
wards 1900 kms -1 at +14° 41' where it is suddenly arrested. 
Then at +14° 46' it is abruptly wrenched upwards again to- 
wards 2100 kms -1 at 14° 49'. 

Fig. 2 is a blow-up of the source region superimposed 
on a far deeper CCD opti cal image, illustratin g that there 
is no optical counterpart jMinchin et al] 1200511 . An optical 
spectrum from the 6.5-m MMT in Arizona of the small, faint 
galaxy 'A', which is superposed on the highest Hi contour 
(1 x 10 20 cm -2 ) at declination +14°47.4', shows that it is 
at a redshift of z = 0.25 and is therefore unconnected with 
VIRGOHI 21. The 17th magnitude galaxy 'C to the left at 
+ 14° 45' is an Hi point-source at this resolution. By compar- 
ison, VIRGOHI 21 is an extended structure in both dimen- 
sions rather than a collection of discrete compact clouds. 
The velocity-declination plot (right) shows the complex 
kinematic structure of VIRGOHI 21. The most remarkable 
feature is the tilted portion between +14°46' (1900 kms -1 ) 
and +14°49' (2100 kms - ) which resembles the signature of 
an ed ge-on rotating disc llKregel. van der Kruit fc de Blokl 
2004), which is what we take it to be. Gas further to the 



North and in particular to the South could either be part 
of the interaction or be clumps further out in the disc. Note 
the lack of connection between galaxy 'C and VIRGOHI 21. 



3 DISCUSSION 

Combining the new data with the old, we argue: 

(a) If attributed to gravitation, changes in velocity of 
galactic size over galactic scales, as seen here, require masses 
of galactic proportions. On dimensional grounds the veloc- 
ity wrench AV (~ 200 kms -1 ) seen at VIRGOHI 21 over a 
conservative length-scale Ax (~ 14 kpc ~ 5 x 10 22 cm ~ 3 ar- 
cmin at the Virgo Cluster distance of 16 Mpc jMinchin et alJ 
120051) 1 implies a mass M > (AV) 2 Ax/G ~ 10 10-11 so- 
lar masses. More specifically the gravitational (free-fall) 
timescale (Ax) 3/2 /(GM) 1/2 needs to be < Ax/V, the time 
it takes for gas travelling at V (relative to M) to change 
velocity by AV. And where AV ~ V, as it appears to 
be here (see Fig. 2b), M > {AV) 2 Ax/G again. For non- 
gravitational alternatives see below. 

(b) The bridge between NGC 4254 and VIRGOHI 21 re- 
veals that they have interacted. Accordin g to published nu- 
merical simulations dVollmer et al.ll2005ft (below) the mor- 
phological peculiarities of NGC 4254 require a perturbing 
mass of ~ lO n M0 - providing independent evidence that 
VIRGOHI 21 weighs ~ lO n M . The simulation implies an 
interaction 3 x 10 s years ago. As the projected length of the 
bridge is 120 kpc this would imply it has been drawn out at 
a projected speed of 390 kms -1 , which is comparable to the 
radial velocity difference between NGC 4254 and VIRGOHI 
21 of 400 kms -1 . 

(c) NGC 4262, the spiral to the North-East, is not in- 
volved. There is no bridge to it (see animation) and the 
radial velocity difference between NGC 4254 and NGC 4262 
is too large (900 kms -1 ) to generate one. 

(d) If object C, the Hi galaxy just to the east (left) of 
VIRGOHI 21 (Fig. 2) were involved its mass, by (a), must be 
~ 1O 11 M0. However its measured Hi velocity distribution 
and size suggest a probable mass < 1O 9 M0, while it has a 
luminosity of only 10 s Lq and a Mhi of 4 x lO 6 Af . It is two 
orders of magnitude too undermassive and underluminous to 
explain the dynamics of VIRGOHI 21, while there is no sign 
of its interacting. 

(e) The structure of VIRGOHI 21 (Fig. 2a) is centred at 
R.A. 12 h 17 m 52 a and elongated in the North-South plane be- 
tween + 14°4l' and +14° 49'. Its remarkable velocity change 
(from 1900 to 2100 kms -1 ) is, however, confined to the tilted 
structure (Fig. 2b) between +14°46' and +14°49', which 
closel y resembles th e characteristic shape of an edge-on disc 
llKregel et al] 12004ft . If indeed what we are witnessing here 
is a dark, gravitationally bound, edge-on rotating disc then 
its properties are as presented in Tabled with a minimum 
mass of 2 x lO lo A/0. Judging from visible disc galaxies, 
whos e masses continue to rise beyond their H I edges dSaluccil 
1997), the full size and mass of such a disc could easily reach 
~ 1O 11 M0. The very low s urface-brightness lim its (dimmer 
than 27.5 B mag arcsec -2 dMinchin et al. 2005)) imply that 
the disc has a Mdynj 'Lb ratio of at least 750 Mq/Lq where 
normal galaxies have < 50. Deep HST observations, capable 
of reaching individual Red Giant stars in VIRGOHI 21, are 
in hand. 
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Figure 1. (a) Hi contour map of the 21-cm observations, superimposed on a 1 square degree 
negative Digitized Sky Survey image. Contours are from 2.5 X 10 19 to 2 X 10 20 cm at intervals 
of 2.5 X 10 19 cm - 2 .(b) Shows the declination-velocity projection of the data cube. More detail 
can be seen in the on-line animation of the cube. 




Table 1. PROPERTIES OF THE DARK DISC 

Diameter 2R, [from +14°46' to 49'] 14 kpc 

Circular Velocity V c , [(2100 - 1900)/2] 100 kms" 1 

Spin Period P, [2nR/V c ] 4 X 10 8 years 

Total Mass M T , [RV^/G] 2 x 1O 1O M 

Face-on Mass-density [M t /ttR 2 ] 2 X 10~ 2 g cm" 

Hydrogen Mass M HI , [F HI = 0.7 Jykms" 1 ] 4 X 1O 7 M 

Face-on gas density Nhi, [Mhi /^rnuR 2 ] 3 X 10 19 cm~ 2 

Total Mass to Blue Light Ratio M HI /L B > 75OM /Lq 
[assuming a disc 0.5' X 3'1 



The integrated spectrum, minus the bridge, is consistent 
with earlier single-dish observations, implying that the in- 
terferometer misses little H I. The surprising ease with which 
it has been mapped is due to two pieces of luck: it is edge 
on, thus increasing the apparent surface density, and most of 
its gas is spread over fairly low velocity widths within each 
synthesised beam. 

The incontrovertible evidence (see animation) of an 
interaction with NGC 4254 provides independent support 
for the massive nature of VIRGOHI 21. NGC 4254 is a 



luminous one-armed spiral gal axy sufficiently peculiar to 
have attracted several studies ^ yejOkamujcLfc Watanabe 
119821: iPhookun. Voeel fc Mundvlll993llVollmer et alJl2005rT 
Single-armed spirals are invariably the result of in teractions 
with close- by massive companions lllve et alJl982l) . The lack 
of any visible companion thus triggered obser vations and dy- 
namic al models. Recent numerical models bv lVollmer et all 
(2005) indicate that NGC 4254 "had a close and rapid en- 
counter with a lO n M0 galaxy ~ 250 Myr ago. The tidal 
intera ction caused the spiral structure...". IPhookun et alJ 
(1993), in a VLA study of the galaxy, find a trail of gas 
leading away from it (their Figure 5) in both the right di- 
rection and with exactly the right velocity-gradient required 
to intersect VIRGOHI 21. Thus the case for VIRGOHI 21 
being the aforesaid ~ 1O 11 M mass which caused the pe- 
culiarities in NGC 4254 seems strong. A detailed simulation 
which includes our new data could remove any doubt. 

Models for VIRGOHI 21 now have a number of crucial 
observations to explain: the broad velocity- width in a galac- 
tic volume (implying large mass); its elongated geometry; 
its steep velocity profile; the bridge to NGC 4254, and to no 
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other galaxy; the damage to NGC 4254 (large mass); and 
the lack of light. Given the above observations we now dis- 
cuss a number of hypotheses as to the origin of VIRGOHI 
21 and its associated Hi. 

(a) We are detecting tidal debris left by the past 
interaction of NGC 4254 and another galaxy. 

This is rather easy to dismiss here where the lines are 
so broad, and there are no such interactors visible in the 
immediate vicinity. Imagine two galaxies with radial veloci- 
ties, Vi and V2 at either end of an approximately linear tidal 
bridge of physical length d pitched at an angle 9 to the plane 
of the sky. A telescope pointed towards it has a transverse 
beam diameter of 6 at the bridge. The only significant gas 
motions within the bridge will be streaming velocities along 
its length. From end to end of the bridge the radial velocity 
difference is \V2 — Vi\ while within the telescope beam the 
measured velocity width, AV, will be (b/dsinO) x | Va — Vi\. 
But, as is well known, bridges of any size arise only when 
the total velocity difference between the interacting galaxies 
- [ V2 — Vi|/cos# here - is of the same order as the ci rcu- 
lar velocity V c in the donor <|Toomre fc Toomrel ll972Tl . It 
follows immediately that AV/V C — (fo/d)tan#. Thus broad 
line widths AV ~ V c , as here, can only be seen within a 
beam if both interactors appear to lie within, or very close to 
the beam (dcos# < b). There are no such putative interac- 
tors even within the Arecibo beam (b ~ 3!6) or close to it, 
otherwise we would see them in Fig. 2a. VIRGOHI 21 cannot 
be such tidal debris. And one cannot escape this conclusion 
by presupposing a very ancient interaction almost along the 
line -of-sight (8 ~ 90°) - if so, wh e re is th e culprit? 

iBekki. Koribalski fc Kilbornl <2005bl) have carried out 
numerical simulations to model interactions that might lead 
to tidal features detectable at 21 cm. They have proposed an 
interaction with NGC 4254 as the likely cause of VIRGOHI 
21. Their Model 1, which does create a cloud of the right 
velocity width - because of its projection along the line of 
sight, fails to show the perturbing galaxy, and in any case 
would never have been claimed as a plausible dark-galaxy 
candidate by us because it would fail the str ingent 'timing- 
argum ent' explained in our previous paper (IMinchin et alJ 
2005). Additionally, the sense of their velocity field relative 
to the donor, NGC 4254, is opposite to the sense actually 
observed. Their favoured Model 4, which does indeed recre- 
ate the disturbance of NGC 4254, has a velocity width of 
20, not 200, kms -1 ! In fact, their simulations demonstrate 
just how hard it is to explain VIRGOHI 21 as tidal debris. 

Other hypotheses have also been put forward, and 
should be examined to see if they can explain our obser- 
vations: 

(b) Two superposed Hi Clouds. 

The components of VIRGOHI 21 are connected both 
spatially and in velocity, making it exceedingly improbably 
that they could be chance superpositions of clouds, while 
the bridge to NGC 4254 is not explained by this hypothesis. 

(c) A tidal tail from a galaxy merging with NGC 
4254; 

iPhookun et al] (I1993T) suggested that the distortion 
of NGC 4254 could be due to infalling gas-clouds. 
Could VIRGOHI 21 be a ta i l, sim ilar to those seen 
in UGC 10214 jBriggs et al] 1200 ll) or NGC 4038/9 
jGordon. Koribalski fc Jones! I 2001]) ? In NGC 4038/9, in 
particular, there is what appears to be a tidal dwarf near 



a strong concentration of H I at the end of the southern tail, 
around 60 k pc from the centre of the system (for a distance 
of 13.8 Mpc: ISaviane. Hibbard fc Richil 2004h . 

Although this hypothesis can, at first sight, explain the 
WSRT observations, with the bridge as an H I tail and VIR- 
GOHI 21 as a tidal dwarf forming in a concentration at the 
tip, it is clear on more detailed examination that this cannot 
be the case. The velocity field in VIRGOHI 21 changes di- 
rection in a way not seen in examples of tidal dwarfs, and to 
a much greater extent - the tidal dwarf in the NGC 4038/9 
system, for instance, has a gradient in the same sense as 
the gas in the tail nearby, whilst VIRGOHI 21's gradient is 
in the opposite direction and has twice the velocity width. 
Even more troubling for this hypothesis, NGC 4254 does not 
show signs of having recently merged with another L* sized 
galaxy - and yet it would have had to do so to have thrown 
out a tidal tail twice as long as that from the violently in- 
teracting NGC 4038/9 system. 

(d) A high column-density part of a giant H I 
rin g. 

Bckk i et a ll (l2005aTl propose that some H I clouds with- 
out optical counterparts could be the high-density regions 
of Hi rings. In their sc enario, objects such as the Leo Ring 
JSchneider et aljll983l) are formed by the tidal stripping of 
low surface-brightness galaxies with extended gas disks. As 
the column-density varies around the ring, it is possible that 
only part of it will be detectable - possibly at a large dis- 
tance from the original galaxy from which the gas has been 
stripped. This would then be identified as an intergalactic 
Hi cloud. 

While this hypothesis might have explained the single- 
dish observations, it is hard to make the WSRT data fit. 
The ring should be orbiting around a large centre of mass 
to the east of VIRGOHI 21 and at a higher velocity, which 
cannot be identified. Also, if the bridge is part of the ring, 
then it is unconnected with the perturbation of NGC 4254 - 
which is both very unlikely and leaves the single-arm mode 
of that galaxy unexplained. 

(e) A three-body interaction. 

It might be that a third mass, possibly object C, has in- 
terfered with the interaction between NGC 42 54 and another 
body (possibly NGC 4262, as proposed by IVollmer et al] 
(2005), as there are no other obvious candidates). However, 
it seems very unlikely that C has sufficient mass to cause 
a perturbation of over 200 kms" 1 to the tidal stream, and 
VIRGOHI 21 passes to its west, not between it and NGC 
4262 as might be expected if it had pulled the stream west- 
ward. The radial velocity of C means that it must be moving 
past the stream at a velocity (relative to the putative undis- 
turbed velocity of the stream at that point of 2100 kms - 
- the velocity of the high-velocity end of VIRGOHI 21) of 
at least 350 kms - . At this speed, it would not have stayed 
close to the stream long enough to have severely disturbed 
it. 

It would also be expected that if C were involved in 
the interaction then there would be gas falling onto it, but 
this does not appear to be the case as there is no gas seen 
between VIRGOHI 21 and C. Nor does C show any signs of 
being disturbed itself; its H I is unresolved, implying that it 
is confined to the area of the optical galaxy, and its optical 
image similarly shows no sign of any perturbation. 

(f) Ram pressure stripping. 
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Uosterloo & van (jorkorn ( 2000 J argue that another ri I 
cloud in the Virgo Cluster, VIRGOHI 4 dPavies et all2004) . 
is caused by ram-pressure stripping from NGC 4388 due to 
an interaction with the hot-gas halo of the M86 sub-group 
and suggest a similar origin could be possible for VIRGOHI 
21. This would explain the bridge without the need to in- 
voke a second galaxy, either interacting or merging. How- 
ever, this cannot give the steep, reversed velocity gradient 
seen in VIRGOHI 21, nor does it explain the distortion to 
the optical disk of NGC 4254. Ram pressure stripping does 
not, therefore, appear to work as an explanation of this sys- 
tem. 

(g) Harassment. 

Harassmen t occurs in much faster interactions than 
tidal stripping 



ijMoore. Lake fc KatzHl998Tl . thus it offers a 
possible mechanism for removing gas from NGC 4254 while 
leaving the interactor sufficiently distant that it can not be 
easily identified. However, such an interaction, while it might 
remove gas, would not form the tidal bridge seen here be- 
tween VIRGOHI 21 and NGC 4254. 

Our conclusion is that the most plausible explanation is 
that VIRGOHI 21 is the cause of the tidal disturbance not 
the effect of other physical processes. 

As might be expected for a dark galaxies, VIRGOHI 21 
has a shortage of baryons compared to visible galaxies. Visi- 
ble disc galaxies typically have, in the form of stars and gas, 
~ 10 per cent of their dynamical mass in baryons, a figure 
consistent with cosmological models which predict ~ 7 times 
as much Dark Matter as baryonic in the Universe at large 
l|Samccilll997h . Had VIRGOHI 21 10 per cent of its mass 
in baryons it would contain ~ 2 x 10 9 solar masses which 
would, if it were once in the form of H I spread evenly across 
the existing dynamical disc, yield a column density ~ 10 21 
atoms cm -2 - sufficie ntly above the "star-format ion thresh- 
old" of ~ 10 20 cm" 2 (iMartin fe Kennicutdl200lh . or above 
the Too mre criterion fo r gravitational i nstability in a rotat- 
mg disc llToomre j ll964D . to form stars lljimenez et alJll997t 
IVerde etedTl2002 ) . Since it apparently contains no stars it 
must have and must always have had, ten times less baryons 
- a p henomenon often apparent in dwarf galaxies jMateol 
I1998T) . Thus 90 per cent of the baryons one might have ex- 
pected to find in a typical galaxy appear to be missing. 
Whatever the reason, the lack of baryons in such massive 
objects can only decrease their detectability. 



4 CONCLUSIONS 

Models for VIRGOHI 21 now have a number of crucial ob- 
servations to explain: the broad velocity- width in a galactic 
volume (implying large mass); its elongated geometry; its 
steep velocity profiles; the bridge to NGC 4254, and to no 
other galaxy; the damage to NGC 4254 (large mass); and 
the lack of light. Of other models that have been suggested, 
or that we have been able to devise, neither tidal debris, nor 
Ram Pressure Stripping nor a Group gravitational field, nor 
long-range interactions deliberately aimed at simulating this 
object can explain the abrupt changes of velocity. 

The new observations make it even harder to escape the 
inference that VIRGOHI 21 contains a massive dark disc. 
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